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ABSTRACT 



Wave propagation in a random, inhomogeneous ocean is 
treated as transmission thru a linear, time-variant, space- 
variant, random communication channel. A consistent notation 
(vis-a-vis ad hoc), fundamental input-output relations, and 
various time-space transformations for both deterministic and 
random linear, time-variant, space-variant, filters are esta- 
blished. Using the method of separation of variables and the 
W.K.B. approximation, a time-invariant, space-variant, random 
transfer function of the ocean volume is derived . The ocean 
volume is characterized by a random index of refraction which 
is a function of depth. The index of refraction is decomposed 
into a deterministic component and a zero mean random compo- 
nent. In addition, two example calculations are made. The 
first example involves the derivation of the equations for 
the random, output electrical signals at each element in a 
receive planar array of complex weighted point sources in 
terms of the frequency spectrum of the transmitted electrical 
signal, the transmit and receive arrays, and the transfer 
function of the ocean medium. The second example involves 
the derivation of the coherence function, i.e., the auto- 



correlation function of the transfer function. 
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I. INTRODUCTION 

Since the wave equation for small amplitude acoustic sig- 
nals is linear , we can represent the ocean medium as a linear , 
time-variant, space-variant, random filter (system or communi- 
cation channel) in general. With this interpretation in mind, 
refer to Fig. 1 which illustrates a basic bistatic communica- 
tion channel geometry, and Fig. 2, which is a mathematical 
block diagram representation of Fig. 1. With respect to Fig. 1, 
both the transmit and receive apertures (arrays) are, in gen- 
eral, volume apertures (arrays) and in motion. Before pro- 
ceeding further, a word of caution concerning Fig. 2. Note, 
for example, that X^ ^ XD T , ^ an< ^ Y ^ Y M°R an 9 enera ^-- 

The equations required for describing the filter's input-output 
relationships and for coupling the transmitted and received 
electrical signals to the medium via the transmit and receive 
apertures are developed in Sections II and III, respectively. 

Let us now describe the notation used in Fig. 2. The posi- 
tion vectors r and r refer to spatial coordinates (x ,y ,z ) 

— o — o o o 

and (x,y,z), respectively, and t refers to time in sec. The 
parameters f and n are frequencies in HZ. where f represents 
input or transmitted frequencies while n represents output or 
received frequencies. Note, that if n ^ f, Doppler spread is 
implied . 

The quantities a, y_, and y_ are vectors whose components 
are spatial frequencies with units of cycles/m. Since spatial 
frequencies are related to both direction cosines and wavelength. 
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Fig. 1. Basic bistatic communication channel 
geometry. 
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Fig. 2. Mathematical block diagram representation of the basic 
bistatic communication channel geometry depicted in 
Fig. 1. 
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and hence, wavenumber components, they represent directions 
of wave propagation. The vector v_ represents input or trans 
mitted spatial frequencies into the medium as in X M (f,v_), 
while represents output or received spatial frequencies 
from the medium as in (n ,S_) . Note, that if jL ^ angular 
spread (scatter) is implied. 

The remaining expressions found in Fig. 2 are further 
described in the following list: 



x ( t , r ) - 


input electrical signal to transmit electro- 
acoustic transducer applied at time t and 
spatial location r of transducer. 


X(f,a) - 


frequency (f) and angular (<*_) spectrum of input 
electrical signal. 


A^, ( f , r ) - 


complex frequency response at spatial location 
r of transmit transducer. Also referred to as 
the complex transmit aperture. 


D T (f,a)- 


transmit far-field directivity function or beam 
pattern. 




input acoustic signal to the medium applied at 
time t and spatial location r. Also, output 
acoustic signal from transmit electro-acoustic 
transducer. 




frequency (f) and angular (v.) spectrum of input 
acoustic signal. 


V T 'V 


t,r) - time-variant, space-variant impulse re- 
sponse of the ocean medium. It represents the 
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response of the medium at time t and spatial 
location r due to the application of an unit 
impulse at time Ct-x)sec,, or t sec. ago, at 
a distance Ir-r^J m. away (see Fig. 1). 

H M (f,_v;t,r) - time-variant, space-variant transfer 
function of the ocean medium. 

y M (t,r) - output acoustic signal from the medium at time 
t and spatial location r. Also, input acoustic 
signal to receive electro-acoustic transducer. 

Y M (n,_S) - frequency (n) and angular (B_) spectrum of out- 
put acoustic signal. 

A^Crift) - complex frequency response at spatial location 
r of receive transducer. Also referred to as 
the complex receive aperture. 

D^(n > 0 ) - receive far-field directivity function or beam 
pattern. 

y(t,r) - output electrical signal from receive electro- 
acoustic transducer at time t and spatial loca- 
tion r of transducer. 

Y(n,x) “ frequency (n) and angular (^) spectrum of output 
electrical signal. 

As was mentioned previously, we can represent the ocean 
medium as a linear, time-variant, space-variant, random filter. 
The term " time- variant" implies motion amongst targets, the 
ocean surface, discrete point scatterers, and the transmit and 
receive apertures (arrays) . Discrete point scatterers in the 
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ocean may include, for example, gas bubbles, fish, and other 
particulate matter. The time-variant property results in 
both Doppler spread and spread in round-trip time delay values. 
If the filter is time- invariant, then no motion is implied. 

As a result, there will be no Doppler spread and no spread in 
round-trip time delay. 

The term "space-variant" implies that the sound speed 
profile (index of refraction) of the ocean is a function of 
position. The space-variant property results in scatter or 
angular spread due to refraction. If the filter is space- 
invariant, then an isospeed medium is implied. As a result, 
there will be no refraction, and hence, no scatter or angular 
spread since the sound rays will be travelling in straight 
lines . 

In addition, since any motion and/or the index of refrac- 
tion can be decomposed into a sum of deterministic (average) 
and random (fluctuating) components, these random components 
can be accounted for via a random filter representation vis-a- 
vis a deterministic filter representation. 

By using a systems theory approach, surface, volume, and/or 
bottom reverberation returns can be modelled as the outputs 
from linear filters. In addition, target returns can also be 
modelled as filter outputs. Furthermore, different transmit 
signals and transmit and receive directivity functions can 
easily be coupled to various models (i.e., transfer functions) 
of the random, inhomogeneous ocean medium in a straightf orward 
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and logical fashion in order to study their effects on target 

l 

detection or parameter estimation using various space-time 
signal processing algorithms. 

The approach of treating the ocean as an isospeed medium, 
and hence, as a linear, time-variant, random communication 
channel is well established [1-19]. This linear, time-varying, 
random system theory approach has also been applied to target 
scattering problems in radar astronomy [20] and to communica- 
tion channels in general [21-23]. However, with respect to 
target models, past research efforts have been devoted mainly 
to the slowly fluctuating point target problem [24-31]. Ef- 
forts to treat more complicated target models were made by 
Kooij [32], Moose [10], and Ziomek and Sibul [19,33]. Kooij 
[32] and Moose [10] both modelled the target as a linear, time- 
invariant, deterministic filter while Ziomek and Sibul [19,33] 
modelled the target as a linear, time- varying , random filter. 

In addition, Ziomek [34] has shown that the form of the gen- 
eralized ambiguity function can be derived by treating the 
scattered acoustic pressure field from a point target (in 
relative motion with respect to a bistatic transmit/receive 
array geometry) as the output of a linear, time-varying, random 
filter . 

Some work has been done in treating the ocean medium as a 
linear, time-variant, space-variant , random filter by Laval 
[9,35] and Laval and Labasque [36]. However, the notation 
used to incorporate the space-variant property is ad hoc, i.e., 
spatial variables are simply included in the arguments of the 
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impulse response and transfer functions, for example, rather 
than haying eyolved from a systematic and consistent notation 
based upon linear, time-varying, space-varying system theory. 

In addition, Laval and Labasque [36] assume functional forms 
for the ocean transfer function instead of deriving them. 
Middleton [37,38] also studied underwater acoustic propagation 
in a random, inhomogeneous ocean, but did not concern himself 
directly with the derivation of random, time-variant, space- 
variant ocean transfer functions. He described the propagation 
phenomena using space- time operators. 

In this paper we will study underwater acoustic propagation 
in a random, inhomogeneous ocean by treating the ocean medium 
as a linear, time-variant, space-variant, random filter. Sec- 
tion II is devoted to a discussion of the fundamentals of 
linear, time-variant, space-variant filters and is based upon 
the generalization of the results contained in Ziomek [39] . 

A consistent notation is developed in a systematic manner. Va- 
rious system functions are introduced and important input-output 
relations and multidimensional (time- space) Fourier transform 
pairs are derived for both deterministic and random filters. 

To the best of the author's knowledge, the expressions presented 
in Section II have not appeared previously in the literature. 

The equations necessary to couple the medium's transfer 
function to the far-field beam patterns of the transmit and 
receive apertures (arrays) and to the frequency spectrum of 
the transmitted signal are discussed in Section III. 
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Section 4.1 is devoted to the main problem of the paper 
which is the derivation of a random ocean transfer function 
incorporating the W.K.B. approximation. In the process of 
the derivation, the index of refraction is decomposed into 
deterministic and random components and it is assumed that 
the medium is "weakly scattering". The transfer function 
derivation was motivated by the work of Clarke [40] on the 
application of the W.K.B. approximation. 

Section 4.2 is devoted to an example calculation of the 
equations for the random, output electrical signals appearing 
at each element in a receive planar array of complex weighted 
point sources in terms of the frequency spectrum of the trans- 
mitted electrical signal, the transmit and receive arrays, and 
the random ocean medium transfer function. 

Finally, in Section 4.3, the autocorrelation function of 
the transfer function, which is also known as the coherence 
function, is calculated. 
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II. FUNDAMENTALS OF LINEAR TIME-VARIANT SPACE-VARIANT FILTERS 
2 . 1 Deterministic Filters 

2.1.1 Impulse response and transfer functions 
A linear, time-variant, space-variant, filter is 
depicted in Fig. 3, where it is characterized by its corres- 

ponding time-varying, space-varying impulse response Mt,^; 
t,r) . The function Mt,^; t,r) describes the response of the 
filter at time t and spatial location r = (x,y,z) due to the 
application of an unit impulse at time (t-x) , or t seconds ago, 
and at a distance |r-r I meters away where r = (x ,y ,z ). 

Note that 

h(x,r ; t,r) = h(t,r; t-t, r-r ) . (2.1-1) 

— o — — — — o 

The relationship between the input signal x(t,r) 
and the output signal y(t,r) is given by 



y ( t,r) 




r-r ) h ( x , r ; 
— — o — o 



t,r) didr 
— — o 

( 2 . 1 - 2 ) 



where it should be noted that both the input and output signals 
are functions of time and space. 



Example 2.1-1 

Different forms of Eq. (2.1-2) can be obtained by making 
the following simplifying assumptions: 
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x ( t, r 




y ( t, r ) 



Fig. 3. Representation of a linear, 
time-variant, space-variant, 
filter . 
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(1) if the linear filter h is time- invariant and space-inya- 
riant, then 

h ( t - [t-x] , r - [r-r^] ) 

Mt,^) (2.1-3) 

and as a result, Eq. (2.1-2) reduces to 



h(x,r o ; 



t,r) = 



y ( t , r ) 




r-r )h(x,r )dxdr 
— — o — O — ' c 



(2.1-4) 



which is a multidimensional convolution integral as would 
be expected. 

(2) if h is time- invariant but space-variant, then 



h ( x , r^ ; t , r ) 



h ( t - [t-x], r; r-r^) 



= h ( x , r ; r-r^) 

= h ( x , r ; r) (2.1-5) 

— o — 



and as a result, Eq. (2.1-2) reduces to 



00 CD 



y(t,r) = 




— ‘X> — oo 



x(t-x, r-r )h(x,r ; r)dxdr . 

— — o — O — — o 



( 2 . 1 - 6 ) 
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(3) if h is tinje-yariant but space-invariant , then 

h(x,r o ; t,r) = h ( t, r - [r-r^] ; t-x) 

= h(t,r ; t-x) 

— o 

= Mx,^; t) (2.1-7) 

and as a result, Eq. (2.1-2) reduces to 



CO CO 



y (t. 



-// 



x(t-x, r-r^ ) h ( X , r^ ; tjdxdr^. 



— oo — oo 



( 2 . 1 - 8 ) 



***** 



Note that if Eq. (2.1-2) is rewritten as 



y (t. 



OO CO 

/ /' 



x(t-a, r-£)h(a,e; t,r)dad£ (2.1-9) 



— OO —CO 



and if the input is an unit impulse applied at time a = (t-x) 
seconds at a distance |c_| = | r-r^ ] meters away, i.e., if 



x(ct,£) = 5 (a - [t-x], £ - [r-r ] ) , (2.1-10) 



then substituting Eq. (2.1-10) into Eq. (2.1-9) yields 
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y (t,r) 



oo oo 




5 (t - a, 



r -c)hCci,c; 
— o — — 



t, r) dctd£ 



= h(T,r o ; t, r) . (2.1-11) 

Analogous to the frequency response or transfer 
function H(f) of linear, time- invariant systems is the time- 
varying, space-varying, frequency response or transfer function 
H(f,y_; t,r) of linear, time-variant, space-variant, systems. 

It is defined as follows: 

H ( f , v ; t,r) = F F {h(t,r ; t,r)} (2.1-12) 

— — x r — o — 



or 



H ( f , v ; t,r) = 



OO oo 

1 7 “ 

— oo — oo 



(x,r ; t ,r) exp (- j 2^f t ) exp ( + j 2 ttv_* r ) dtdr 



(2.1-13) 



where f corresponds to input frequencies in HZ. and is a 
vector whose components are input spatial frequencies with units 
of cycles/meter. As was previously mentioned in Section 1, 
spatial frequencies are related to wavenumber components, and 
hence, they represent directions of wave propagation. Similarly, 
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h (t , r ; t, r) 
— o — 



F f 1 F~ 1 { H (.f , v_; t, r ) } 



(2.1-14) 



or 

hfx,^; t,r) 

The choice of a plus (+) sign in the exponent of 
exp ( + j 2tiv_* appearing in Eq. (2.1-13), which corresponds to 

the forward spatial Fourier transform w.r.t. r^, was not 
arbitrary. This choice of sign convention is meant to be con- 
sistent with that of the spatial Fourier transform relation- 
ship between a complex aperture function and its directivity 
function (beam pattern) as is developed later in Section 3. 
Besides, the integrand term 

exp (- j2irf x ) exp ( + j27iy r^) 

appearing in Eq . (2.1-13) has the nice physical interpretation 

of being a time-harmonic plane wave travelling in the direction 
of increasing r Q = | J with the sign convention as given. 

Example 2.1-2 

If the linear filter h is time- invariant and space-in- 
variant, then its corresponding transfer function can be ob- 



OO OQ 

// 

-oo —co 



H (f t , r) exp ( + j 2rr fx ) exp (- j 2 ttv_* r^) dfdv_ 



(2.1-15) 
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tained by substituting Eq. (2.1-3) into Eq. (2.1-13). Doing 
so yields the following result; 



H ( f , v ) = 



oo oa 

If 

— oo — oo 



h (t , r^) exp (- j2iTf t ) exp ( + j27rv_* r^) dtdr^ 



(2.1-16) 



Similarly, from Eq. (2.1-15), 



h(T,r o ) 



OO 00 

■II 

— ■ oo — oo 



H (f ,y_) exp (+j 2if x ) exp (- j 2uy_ 



r )dfdv, 
— o - 



(2.1-17) 



* ★ ★ * ★ 



Calculating the transfer function according to 
Eq. (2.1-13) assumes that one knows the impulse response func- 
tion. However, even if the impulse response is not known, 
the transfer function can still be obtained by using a time- 
harmonic plane wave as an input signal. This can be shown by 
representing a linear, time-varying, space-varying, filter by 
the linear operator L(*) which operates on input signals that 
are functions of both time and space. The output of the 
filter y(t,r) can then be expressed as 



y(t,r) = L 



[x (t,r) J . 



(2.1-18) 
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If we let 



x(t,r) = exp (+ j 2tt f t) exp (- j 2ttv_* r) , (2.1-19) 

which is in the form of a time-harmonic plane wave travelling 
in the direction of increasing r = |r|, then from Eqs. (2.1-2) 
and (2.1-18) we can write 



L [exp ( + j 2uf t) exp (- j 2iry_* r) ] 



oo 



— CD 



exp[ + j2^f (t-r)] exp[-j2uv_* (r-r^) ] • 



h (x / r_ ; t,r)didr , 
— o — — o 



( 2 . 1 - 20 ) 



and by using the definition of the transfer function given by 
Eq. (2.1-13), Eq. (2.1-20) can be expressed as 

L [exp ( + j 2 ^f t) exp (- j 2 ttv_* r ) ] = H(f,v_; t , r) exp ( + j 2:rf t) exp (- j 2™* r ) 

( 2 . 1 - 21 ) 

where L [exp ( + j 2^f t) exp (- j 2™ • r ) J is the response of the filter 
to exp ( + j 2iTf t) exp (- j 2 tt v_* r ) , and H(f,v_; t,r) is the time-varying, 
space-varying transfer function of the filter evaluated at 
f and v_. 

Equation (2.1-21) is a fundamental result that 
will be used in the derivation of an ocean transfer function 
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in Section 4.1. In addition, Eq. (.2.1-21) can be used to 
represent the output y(t,r) in terms of the transfer function. 
For example, if we express the input x(t,r) in terms of the 
following multidimensional inverse Fourier transform 



x(t,r) = Fj-V^fXCf , v) } 



( 2 . 1 - 22 ) 



or 



GO CO 



X (t. 



“■// 



X(f , v) exp (+j2irf t) exp (- j 2ity_*r) dfdv_, 

(2.1-23) 



then from Eqs. (2.1-18) and(2.1-21), we can write the output 
y(t,r) as 



y(t,r) = 



00 00 

// 



X (f , v) L [exp (+ j 2rrf t) exp (- j 2frv_*r) ] dfd_v . 

(2. 1-24] 



or 



Y (t,r) 



03 JO 




t, r) exp ( + j 2iTf t) exp (- j 2iry > *r) dfd 



(2.1-25) 



Compare Eq. (2.1-25) with Eq. (2.1-2). 
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2.1.2 Spreading and bi-frequency functions 
Two additional filter functions will now be intro- 
duced; namely, the spreading function and the bi-frequency 
function . The spreading function will be discussed first. 

The spreading function S(x,r <3 ? 4>,>0 is defined as 

S ( t ,r ; <J>,k) = F.F fh(t,r ; t,r)} (2.1-26) 

— o — t r — o — 



or 



S ( x , ; ij) , _< ) = 



oo oo 

J j h(T,r Q ; t , r ) exp (- j 2-rr 4> t ) exp ( + j 2u£* r) dtdr 



— CO -CO 



(2.1-27) 

where $ corresponds to the rate of change of the filter's im- 
pulse response in HZ. and <_ is a vector whose components are 
spatial frequencies which correspond to the rate of change of 
the impulse response in cycles/m. Similarly, 



h (t ,r ; t , r ) = F~ 1 F~ 1 {S (t , r ; <{.,<)} ( 2 . 1 - 28 ) 

— o — <p < — o — 



or 



h ( t , r ; t,r) 
— o — 



oo at 

// 

__ JO — OO 



S(t 






4> , k ) exp ( + j 2 tmj> t) exp (- j 2 tt r) d$d<_. 



(2.1-29) 
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An alternate representation of the output y(t,r) 
can be obtained by substituting Eq. (2.1-29) into Eq. (2.1-2). 
Doing so yields 



y(t,r) = 



oo oo oo oo 

//// 

— oo — oo — oo — oo 



x(t-t, r-r^) exp ( + j 2ir<f>t ) exp (- j 2™*r ) • 

S(t r ; $ , k) d<j)dxdTdr (2.1-30) 

— o — — — o 



where the integrand term 



S(T / r_ o ; <i> , O x ( t-t , r-r^) exp ( + j2ir<|>t) exp (- j2'tf£ , r) 

is a time, space, and "frequency" (frequency in HZ. and spatial 
frequencies in cycles/m.) shifted replica of the input signal 
x(t,r) weighted by the spreading function S(T > r Q ; <!>,£). 
Therefore, for a given input (source) location as specified by 
r^, the spreading function determines the amount of spread in 
round-trip time delay r (sec.) , frequency <j>(HZ.), and spatial 
frequencies £ (cycles/m.) that an input signal will undergo as 
it passes through a linear, time-varying, space-varying, channel. 
Equations (2.1-30) and (2.1-2) both indicate that the output 
y(t,r) at some receiver location r is dependent upon the loca- 
tion of the input (source) via the spreading function 
S(T,r Q ; 4>,£) w.r.t. Eq. (2.1-30), and via the impulse response 
h(x,r ; t,r) w.r.t. Eq. (2.1-2). For example, the existence and 
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extent of SOFAR transmission depends on both the depths of the 
source and receiver [41] . 

The last filter function to be discussed is the 
bi-frequency function. The bi-frequency function B(f,y_; 4>,jc) 
is defined as 



B ( f , v ; 



) 4 



F t F r {H(f ,v; 



t , r ) 



(2.1-31) 



or 



CO CO 



B ( f , = J j H(f,y_; t , r ) exp (- j 2n $ t ) exp ( + j 2tt jc • r ) d tdr . 



•CO — CO 



(2.1-32) 



Similarly, 



H ( f , v ; t,r) = F" 1 F" 1 (B(f,v; , < ) > 

— — $ JC — — 



(2.1-33) 



or 



CO oo 



// 

— oo — oo 



H(f,v_; t,r) = / / 3(f,v_; <)>,£) exp ( + j 2u tf> t ) exp ( -j 2tt<_* r) d<()d<_. 

(2.1-34) 



The term bi-frequency function was originally used to denote 
the appearance of the two frequency variables f and 4> in linear, 
time-varying systems theory. We can generalize this notion 
by using the term bi- frequency function to denote the appearance 
of the input "frequency" pair f and v_, and the filter variation 
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"frequency" pair <j> and <_ in linear, time-yarying , space-yary- 
ing, systems theory. 

Just as the spreading function gives an indication 
of how rapidly the impulse response changes with time and 
space, the bi-frequency function gives an indication of how 
rapidly the transfer function changes with time and space. 



taking the Fourier transform of the spreading function w.r.t. 
t and r^ , i . e . , 



The bi-frequency function can also be obtained by 



B (f , v_; <j>,j<) = F^F r {S(x,r o ; $,£_)} 



( 2 . 1 - 35 ) 



or 



oo 



ao 




j> , k) exp (-j2irfx) exp ( + j 2Try_‘r Q ) dxdr 



( 2 . 1 - 36 ) 



and, similarly 




( 2 . 1 - 37 ) 



or 



oo 



oo 




$ , <_) exp ( + j 2 irf x ) exp (-j 2 it v_*r^) dfdy_ . 



( 2 . 1 - 38 ) 
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In addition, the impulse response and bi-frequency 
functions form a Fourier transform pair, i.e.. 



B(f,v; <J> ,jc) = F T F r F t F r (h(T,r o ; t,r) } (2.1-39) 

— o — 



or 



oo ao oo oo 



B ( f , v ; <p , <_) = 




h(x,r ; t, r ) exp (- j 2irf x) exo ( + j 2 ttv • r ) 

— O — — — o 



— oo — oo — oo — oo 



exp (-j 2ir(f>t) exp ( + j2^£*r) dxd^dtdr 



(2.1-40) 



and, similarly, 



h(x,r o ; t,r) = F^ 1 f” 1 f" 1 F^ 1 { B ( f , v ; ♦,«<)} 



(2.1-41) 



or 



h ( x , ; t , r ) 




,jc) exp ( + j 2wf x) exp (- j2Tryy r ) 



exp ( + j2in{>t) exp (- j 2 ik *r) dfdyd4>di<_. 



(2.1-42) 



The interdependence which exists amongst the four 
filter functions is illustrated in Fig. 4. Any one of 

these functions may be used to define completely a linear, 
time-varying, space-varying, communication channel. 
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SPREADING 

FUNCTION 



TIME- VARYING SPACE- VARYING 
IMPULSE RESPONSE 




TIME-VARYING 
SPACE -VARYING 
TRANSFER 
FUNCTION 



BI-FREOUENCY 

FUNCTION 



Interdependence amongst the four filter 
functions that characterize linear, time- 
varying , space-varying , communication 
channels . 



Fig, 4, 
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2.1.3 Output frequency and angular spectrum 
The output frequency (q ) and angular (B_) spectrum 
Y(q,0) is defined as 



Y(n ,£) = F t F r (y (t,r) ) 



(2.1-43) 



or 



Y(n, 



OO QC 

“ ! // 



y ( t , r) exp (- j 2irn t) exp ( + j 2ir6_* r) dtdr (2.1-44) 



.00 _ CO 



where n corresponds to output frequencies in HZ . and £ is a 
vector whose components are output spatial frequencies with 
units of cycles/meter. Substituting Eq. (2.1-25) into Eq. 
(2. 1-44 ) yields 



Y(n,S) 



CD OO 




_ 00 _ OO 



OO OO 




— OO — OO 



t , r ) exp [- j 2 ir (n-f ) t] * 



exp [+ j 2u ( 8_-v ) *r]dtdrdfdv 



(2.1-45) 



where, from Eq. (2.1-32), it can be seen that the inner multi- 
dimensional integral is equal to B(f,v_; n-f, £-v_) so that 



Y (n , 8 ) 



OO OO 




n-f, 8-\>)dfdv 



— OO — OO 



(2.1-46) 
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which is in the form of a multidimensional "frequency" domain 
convolution integral. Note that <j> = n-f or n ^ f+4>, and that 
jc = or | = 2L + JS* Thus, the output frequencies n in HZ. are 

equal to the sum of the input frequencies f and the variations 
in frequency, <j> , due to the time-varying property of the filter. 
Similarly, the output spatial frequencies 3_ (directions of 
wave propagation) in cycles/m. are equal to the sum of the in- 
put spatial frequencies v_ (directions of propagation of trans- 
mitted plane waves) and the variations in spatial frequency 
(directions of wave propagation) , jc, due to the space-varying 
property of the filter. Equation (2.1-46) demonstrates that a 
linear, time-variant, space-variant, filter will spread the 
input frequency and angular spectrum X(f,y_) in both frequencies 
in HZ. and spatial frequencies in cycles/m. 

Example 2.1-3 

If- the linear filter h is time- invariant and space-in- 
variant, then H(f,y; t,r) = H(f,y) [see Eq. (2.1-16)] and, as a 
result, Eq. (2.1-32) reduces to 
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B(f,y; $,k) = H(f,v) 



/ 



exp (- j 2 tt$ t) dt 



— oo 



— oo 



or 



B ( f , v ; $, k) = H (f , v) 6 U) 5 (<) 



(2.1-47) 
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since 

F t {l> = 5 (<j>) 

and 

F r {l} = 5(<) . 



Substituting Eq. (2.1-47) into Eq. (2.1-46) yields 



Y(n,B) 



oo no 




X ( f , v ) H ( f , ' j ) 6 ( n — f ) <5 (B-v)dfdv 



which simplifies to 

Y(n,£) = X(n,£)H(n, £) , 

and by replacing n with f and £ with \^, we finally obtain 

Y(f,_y) = X(f ,v)H (f r±) (2.1-48) 

which is the output frequency and angular spectrum from a linear, 
time- invariant , space- invariant , filter. Note that the output 
frequency f and the output spatial frequency vector are iden- 
tical with those of the input X(f,_v). Hence, as would be ex- 
pected, there is no "frequency" spreading. This would corres- 
pond to the physical situation of transmitting a signal via a 
transmit aperture (array) to a receive aperture (array) when 
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the platforms containing the apertures (arrays) are not in 
motion, and the intervening ocean medium has a constant speed 
of sound (index of refraction) and no discrete point scatterers 
in motion. No motion implies no Doppler (frequency) spread 
and a constant speed of sound implies that the sound rays 
travel in straight lines, i.e., there is no angular spread 
(scatter) or change in the direction of propagation of the 
transmitted sound rays. 

★ * ★ ★ * 
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2 . 2 Random Filters 

2.2.1 Filter autocorrelation functions 
In Section 2.1, four system functions were intro- 
duced which are used to characterize linear, time-variant, 
space-variant, filters. However, if the filter is random, 
then each of these system functions must be considered as 
random functions. As a result, we must work with the respec- 
tive system autocorrelation functions which are defined as 
follows : 

R h ( T ' T ' '^-o' ^o : * EthU,^; t, r) h* ( t ' , r^; t',r')} 

( 2 . 2 - 1 ) 

R H (f,f*,v,v'; t, t ' , r , r ' ) = E{H (f , v; t,r) H* (f ' , v ' ; t',r')> 

( 2 . 2 - 2 ) 



Rq ( t , t , r _ , r t $ , $ , k , < ) — E{S(x, r ^ t 1 / $ , ^ ) / 



— o — o 



— o 



(2.2-3) 



and 



R B (f ,f ’ ,v, y' ; <f>,<f>* ,k,k ' ) = E{B(f,y; $ ,k)B* ( f ' ,y' ; 

(2.2-4) 



where E{ • } is the expectation or ensemble average operator 
and the asterisk denotes complex conjugate. 

If we use the sign convention that forward trans- 
forms w.r.t. x,t,r^, and r' are defined with a minus sign in 
the exponent of the complex exponential ( inverse transforms 
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w.r.t. f,<j>,v_', and ic* are defined with a plus sign) and for- 
ward transforms w.r.t. t',t',r , and r are defined with a 
— o — 

plus sign ( inverse transforms w.r.t. f',4>',£, and <_ are defined 
with a minus sign) , then it can be shown that the four system 
autocorrelation functions are related by the following Fourier 
transform pairs: 



V f , 



V , V 



t,t' , r , r 1 ) a F F 



T T 






<v 



T , T 



'^o'£o : 



t,t- 



(2.2-5) 



and 

\ ( T ' T ' >E 0 'E^' t,t',r,r') 




{R H (f,f 



v,£'; t, t ' , r ,r ' ) 
( 2 . 2 - 6 ) 



(f , f ' ; <j> , ) 



F t F t' F r F r' {R H (f ' f ' 



, v, V ' ; 



t,t' , r , r ' ) } 
(2.2-7) 



and 

Rg (f ,f * , v ' ; t, t' ,r,r ' ) 



’T 1 f~]‘f~ 1 f”| {Rq (f , f 1 



V , V 



<P / $ t £/•£ ) ^ / 
( 2 . 2 - 8 ) 



R Q ( t , t ' , r , r^ ; <j> , <j> ' , k , < ' ) 
and 

R u ( t , t ' , r , r! ; t , t ' , r , r ' ) 
n — o — o — — 



t t ' r r 



{ R h ( f , T 






t, t ’ ,r ,r ’ ) } 



(2.2-9) 




{ Rg ( T , T 



f / 1-1 / $ / $ 1 # £ / £* 

( 2 . 2 - 10 ) 



R B (f ,f ’ ,v,v’ ; 1>’ ,£,<_’) = F^F^.F 



r F r'*- a s (x,T ' '^o'--o ! * ' * 

— o — o 



( 2 . 2 - 11 ) 
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and 

***'»£*!’> = F^F^F^F^fR^f ,f ' ,v,v* ; 

( 2 . 2 - 12 ) 



and finally. 



Rg(f ,f ' , v_, v_' ; <P, <P ) = 



F , F ,' F r F r' F t F f F r F r’ (R h (t ' r, '£o'^ : 
— o — o — — 

t , t ’ r , r ’ ) } (2.2-13) 



and 



Bhtt.f.ro.r^, t,f,r,r/) = (f , f ' , v , v • 



^ F 4* f ^ f ^ )}. 



(2.2-14) 



The interdependence which exists amongst the four 
filter autocorrelation functions is illustrated in Fig. 5. 
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2.2.2 Uncorrelated spreading and the scattering 
function 

In this section we will examine the consequences 
of assuming that the spreading function Stx,^; <f>,£) is un - 
correlated with S(t', r^; $',£') for all values of t' / x, 

r ' t* r , <f> ' $ , and < ' ^ < . 

— o — o — — 

The assumption of uncorrelated spreading is mathe 

matically equivalent to stating that the autocovariance of 

S ( x <J>,jc) and S(x',r^; 4 >',£') is zero for all values of 

t’ ? x, r! jt r , <}>' 7 * <}>, and <' / k, i.e., 

— o — o — — 

C_ ( T f T f r / r 7 $ / $ / < / K 1 ) = R q ( X / T / IT / IT 7 $ / $ * / K , 1 C ) “* 

S — O — O — — * 5 — O O “ — 

u s (T,r^; <j),£) <p ' , ) = 0 

x ' t* x , r! r , <ji' 7 * $ , and k ' ^ 
— o — o — 

(2.2-15) 

where Cg is the autocovariance function and $,<_) = 

E(S(x,r ; If it is assumed that y c (x,r ; <}>,<) = 0, 

then Eq. (2.2-15) is equivalent to 

R s ( T / T f r ^ t r^ 7 f $ ) 

where 

2 

R<~» ( x r r / <j) / < ) = E { | S ( x , r 7 <p t k ) I } 



= R 3 ( x , ^7 $ , jc) <5 ( x-x 1 ) * 

5 (r -r • ) 0 (<j>-<j>' ) 6 (k-k: ' ) 

— o — o — — 



(2.2-16) 



( 2 . 2 - 17 ) 
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is called the scattering function and is equal to the mean 
squared value of the spreading function. The scattering func- 
tion can be thought of as an average power density function 
which determines the average amount of spread that an input 
signal's power will undergo as a function of round-trip time 
delay x , frequency <j> , and spatial frequencies <_ for a given 
input (source) location r^. Note that the scattering function 
is a real, positive valued function. 

Equation (2.2-16) is the result of the assumption 
that the spreading function is zero mean. However, if the 
spreading function is non-zero mean, it is convenient to do 
the analysis with the centered process 

S^(T,r ; <fr,K) = S ( T , r ; $,<) - U c (T,r ? *,«). 

C — O — — O — b —O — 

The random process ,r_ Q ; 4>,<_) has zero mean, and as a re- 

sult, its autocovariance function is equal to 



U , T ' 

— o — o 



<j>,4> ' ,K,K ' ) = E( S_ (t ,r ; <j> ,< ) S* (T ' ,£■ ; <t> ' , < ' ) > 

— — C — ' o — c — o — 

s c ~° - 

= Cg ( f , x ' , r^ , r^ ; 4> , 4> *,<,<’) , 



i.e., the autocorrelation function of the zero mean centered 
process is equal to the autocovariance function of the original 
non- zero mean spreading function. 



35 



Let us next examine the effect that the assumption 
of uncorrelated spreading has on the remaining three system 
autocorrelation functions. As Fig. 5 indicates, Rjj can 

be obtained from R g by performing the following transformations: 

Rn^fV^v’ ;t,t'# £,£')= F“ 1 F"lF" 1 F^F x F T ,F r F p , {R s (t,t\ 

? ” — — — o — o 

r o ,r^; * } (2.2-18) 

or 

Rjj (f , f ’ ; t,t’,r,r’) 

exp [+ j 2 .tt ( (j> t-4» ' t ')]exp [- j 2 u (<_•£-<_' ♦£')]• 

exp [ — j 2ir (fx“f'T ' ) ] exp [+ j 2 tt ( r^-v 1 • r^) ] • 

d<j>d(j) ’ d<d< ' dtdt ' dr dr’ . (2.2-19) 

— — — o — o 

If Eq. (2.2-16) is substituted into Eq. (2.2-19), then 

Ry (f , f ' , y_,y’ ; t,t',r,r') = R^AfjAw At,Ar) (2.2-20) 



OO CD 




— OO — OO 



where 



Rjj ( Af , A£; At,Ar) 



CO OO 



— OO — CO 



CO 




OO 



— CD 



R„(-c,r ; <}> , <_) exp ( + j2TT<{>At) • 

exp (- j 2 ttk_* Ar) exp (- j 2iAf x ) • 
exp ( + j 2w A r ) d<{)di<_dTdr , 



( 2 . 2 - 21 ) 
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Af = f-f ' , A_v = v_-y_’ , At = t- 1 ' , and Ar = r-_r ' . 

It can be seen from Eq. (2.2-20) that when uncorrelated spread 
ing is assumed, the autocorrelation function R„ becomes a func 
tion of the differences Af, A_v, At, and Ar only. This implies 
that the random process H(f,_v; t,r) is wide-sense stationary 
in frequency, spatial frequencies, time, and space. An addi- 
tional requirement for H(f,u_; t,r) to be wide-sense stationary 
is that the mean value 

u H (f,£ ; t,r) = E{H(f,v_; t,r) } 

is a constant. Since the four filter functions are related 
by linear transformations (see Fig. 4) , and since it was 

assumed that iig(x ,r^; <}>,<) = 0, then p^tx,^; t,r) = 0, 

= 0, and Pg(f,£; <P , <_) = 0 which are constant, 
zero mean values. Therefore, the condition of uncorrelated 
spreading in round-trip time delay x, input (source) location 
r^, frequency <f> , and spatial frequencies £ is equivalent to a 
condition of wide-sense stationarity in frequency Af, spatial 
frequencies time At, and output (receiver) location Ar, 

respectively. If uncorrelated spreading in x, r , $ , and <_ 

occur together, then we have a wide-sense stationary uncorre - 
lated spreading (WSSUS) communication channel. 

Consider the autocorrelation function R^ next. 

If Eq. (2.2-16) is substituted into Eq. (2.2-10) and the in- 
dicated transformations are performed, then 
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R, (x,T',r , r* ; t,t‘,r,r') = R. (x,r ; At , Ar) 5 (t -t ' ) o (r -r ' ) 
n — o — o — — n — o — — o — o 

( 2 . 2 - 22 ) 

where 



OO CO 



R^ ( t , ; At,Ar) ~ j J R S^ T, ^o ; > *) ex P (+3 2tt <j> a t) • 



— 03 — OO 



exp (- j 2 ti<_* Ar) d4id<_. 



(2.2-23) 



Equation (2.2-22) indicates that the random process h(x,r Q ; t,r) 
is wide-sense stationary in time and space because of the At 
and Ar dependence and since ^(wr^; t/r) = Equation (2.2-22) 

also indicates that t,r) is uncorrelated for all values 

of t ' 7* t ana r^ ^ r^ . 

Finally, if Eq. (2.2-16) is substituted into 
Eq. (2.2-11) and the indicated transformations are performed, 
then 



Rg (f , f ' , y_,v_' ; , <j> ' ,jc,jc ' ) = Rg ( Af , Aw <j> , jc) 6 ( <{>-<{> * ) 6 ( ) 



(2.2-24) 



where 



Rg ( Af , Aw $ , < ) 




R c ( x , r ; <p, <) exp (- j 2irAf x) * 

o O 

exp (+j 2 ttAv • r )dtdr . 
r ■“ — — o — o 



(2.2-25) 



Equation (2.2-24) indicates that the random process B(f,w <f>,<_) 
is wide-sense stationary in frequency and spatial frequencies 
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because of the Af and A\^ dependence and since M B (f,v_; = 0 

Equation (2.2-24) also indicates that B(f,_v; <)>,£) is uncorre- 
lated for all values of <f>’ ^ <P and je' ^ k. 

Therefore, in suirunary, under the assumption of 
uncorrelated spreading the four filter autocorrelation func- 
tions reduce as follows: 



R s < T / T ' /£o'^o ; ) = V T, £o ; 5 ( T_T ' ) 5 






(2.2-16) 



R^f ,f ' ,\>,v/ ; t,t' ,r,r' ) = RjjCAf^v^; At, Ar) (2.2-20) 



*h 



(T,T',r ,£•; t, t ' , r , r ' ) 
— o — o — — 



R.(x,r ; At, Ar) 6 (t-t ' ) 6 (r -r! ) 
n — o — — o — o 

( 2 . 2 - 22 ) 



and 



R B (f ,f • ,v,v* ; *,*',£,£') = R^Af ,Av; 4 >,k) 6 (*-*•) 5 (<-<’) 

(2.2-24) 



where R c (x,r ; <&,<) is the scattering function. 

b “O ” 
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2.2.3 The scattering function and its Fourier 
transforms 

By inspecting Eq. (2.2-23), we can write that 



and 

or 

R s ^ T ' £o 



R h ( t , r ; At, Ar) 



F" 1 F" 1 {R„ (t ,r ; <(.,<) } 

(D < O —O — 



(2.2-26) 



^ s ( T,ro ; *,<) = F At F Ar {R h (T '£o ; At ' A * )} 



(2.2-27) 



CD CD 



9 <) = 



// v " fc 



; At, Ar) exp (- j 2 tti{) At) 



— CO *-oo 



exp (+j 2 ttic • Ar) dAtdAr . 



(2.2-28) 



that 



and 



or 






Also, by inspecting Eq. (2.2-25), we can write 

Rq ( Af , Ay ; <(> ,*_) = R x R r . { R s ^ T ' £o ' ^ } 

— o 

(2.2-29) 

R s (t,r o ; ♦,<) = F"^F~J{R B (Af , Av ; *,<)} 

(2.2-30) 



9 , <) 



CD QQ 




R ( A f , A v ; 4> , ic ) exp ( + j 2 v A f t ) • 

B — — 

exp ( - j 2 tt Ay r )dAfdAy. 



(2.2-31) 
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Additional transform pairs can be 
follows. With the use of Eq. (2.2-23), Eq. (2 
expressed as 



Rjj(Af,A_v ; At,Ar) = F^F r 



or 



CO oo 



Rjj ( Af , iv ; At,Ar) = 



j I R h (t '£o ; 



At, A£ ) exp ( 



— oo — oo 



exp ( + j 2ir Av • r ) drdr 
— — o — o 



and 



R h (T '^0 '* At ' A ^ } = F If F Ai {R H (Af ' A ^ 



or 



oo oo 



R*( t , r ; At, Ar) = 
n — o “ 



r) = j j R^AffAv ; At,Ar)exp 
exp ( - j 2 ;t A£* r )dAfdAv 



— oo — oo 



Next, if Eq. (2.2-28) is substituted into Eq. 



RgUffAjy ; 4>,<) - F t F r F At F A r {R h {t 



obtained as 
.2-21) can be 

; At, Ar ) } 

(2.2-32) 

- j2irAf t)* 

(2.2-33) 

; At, Ar ) } 

(2.2-34) 

( + j 2irAf x ) • 

(2.2-35) 

(2.2-25), then 

At, Ar) } 



(2.2-36) 
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or 



V 4£ , Av ; $,£) 



OO QO OO OO 




; At , Ar ) exp (- j 2:rAf t ) • 

exp ( + j 2 ttAv/ r^ exp (- j 27nj>At) • 
exp ( + j2ir<_* Ar) dtdr^dAtdAr 



(2.2-37) 



and 

R h (t,r o ; At, Ar) = ^F^F^F^ Rg (Af , Av ; <}>,<)} 

(2.2-38) 



or 



V T 



'Ad 



At, Ar ) 



OO OO OO QO 

V 4f , Av ; 4> , O exp ( + j 2irAf t ) • 

exp (- j 2irAyy r^) exp ( + j 2 it 4> At) • 
exp (- j 2 ttjc • Ar ) dAf dAv_d<pdic . 




(2.2-39) 



Finally, with the use of Eq . (2.2-33), Eq. (2.2-37) can be 

expressed as 



Rb< 4£ ,A_v ; <J>,0 = F Afc F Ar {R h ( Af , Av ; At,Ar)} 

(2.2-40) 

or 



Rq ( Af , Av_ ; <J>,jc) 



— CO 



— OO 



R ( Af , A v ; At , Ar ) exp (- j 2 tt <(> A t ) • 
hi — 

exp ( + j 2n<_’ Ar ) dAtdAr (2.2-41) 
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and 

Rjj (Af , Av_ ; At,Ar) = F^F^ ^ { Rg ( Af , Ay_ ; <J> , k_) } 

(2.2-42) 

or 



Rjj(Af,Av ; At , Ar) 



CO CO 




Rg(Af,Av_ ; (j) , ic ) exp ( + j 2ir<j>At) • 
exp (-j2i:K_* Ar) di}idK_. (2.2-43) 



The various Fourier transform pairs are summarized 
in Fig. 6. The scattering function, or any of its Fourier 

transforms, is a complete second order statistical description 
of a WSSUS communication channel. 

Another very important Fourier transform relation- 
ship can be obtained by substituting Eq. (2.2-35) into Eq. 
(2.2-28) which yields 



R S (, '£o 



4> > £.) 



F F 

At Ar~ Af Av 




( Af , A v 



At, Ar) } 

(2.2-44) 



or 



R^ ( t , r Q , (f> , k ) 



00 00 00 oo 

R„ ( A f , A v ; At , Ar ) exp (- j 2 -rr A t ) • 

n — — 

exp ( + j 2 itk_* Ar ) exp ( + j 2-xAf x) • 
exp (- j 2 it A yr^) dAtdArdAfdAv_ 




(2.2-45) 



and 



At, A r ) 
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Fig. 6. The scattering function 
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Rjj ( Af , Av_ ; At,Ar) = 



{R S (, '^o ; *-i )} 

— — o 



( 2 . 2 - 46 ) 



sxnce 



oo oo oo oo 



R H (Af,Av_ ; At,Ar) = 




*.00 ib CD .CO .CO 



Rg(T ,£ o ; p , £_) exp ( + j 2Ti(j)At) ♦ 

exp (- j2ir<_* Ar) exp (-j 2nAf t ) • 
exp ( + j 2 tt Av_* r^) dpd^drdr ^ . 

( 2 . 2 - 21 ) 
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2.2.4 Input-output relations 

We will now proceed to derive an expression for 
the autocorrelation function of the output from a linear, 
time-varying, space-varying, random filter. The output auto- 
correlation function (t, t' ,r,r ' ) is defined as follows: 



yt,t' 


,£,£’) 4 E{y(t',£)y*(t £') }. (2.2-47) 


If Eq. (2.1-30) is 


substituted into Eq. (2.2-47), then the 


output correlation 


function in time and space is given by 



oo 

/° 


00 


Ry ( t , t ' ,£,£') = / 

J 


••• / x ( t— t , r-r ) x* (t ’ - t ' , r'-rl)* 

I — — o — ™0 

J 


— OO 


— oo 

ex P [■*• j 2 TT (4>t-$ ' t')]exp [ — j 2 7T (_<•£ “ £'•£')] 

R S , T ' ,£o' ; ) d 4> d £.dTd£o* 

d(j> ’ die ’ dr ’ dr! . (2.2-48) 

— — o 



If the random filter is a WSSUS (wide-sense stationary uncor- 
related spreading) communication channel, then the autocorre- 
lation function of the spreading function is given by Eq. 
(2.2-16). Substituting Eq. (2.2-16) into Eq. (2.2-43) yields 
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R y (.t, t' ,r,r' ) 



30 OO OO 30 




x Ct-T , r-r ) x*(t'-t , r ' -r )• 
o — — o 



exp ( + j 2 TTij) At) exp (- j 2 tt£ * Ar ). • 
R_(t,r ; <j>, <)d<j>d<dTdr . 

b O O 



(2.2-49) 



The mean squared value of y(t,r) is given by 
E( |y (t, r) | 2 } = R y (t,t,r,r) 
where, from Eq. (2.2-48), we have in general that 



R y (t, t,r, r) 




oo 

/ 



x(t-T,r-r ) x* ( t-t ' , r-r ' ) 
— — o — — o 



exp ( + j 2irA<j>t) exp (- j 2:r A<yr ) • 

R c (T,x',r ,r' ; , <{> 1 ,<,<') d^dtcdtdr • 

^ — o o o 

d$ ' d< ' dt ’ dr ' (2.2-50) 

— — o 



where A4» = and Ai< = <-<_'• In the case of a WSSUS commu- 

nication channel, we have from Eq. (2.2-49) that 
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0O CD OO 0O 



(t,t,r,r) = J j j j |x 



( t- t , r- ) | Rg ( t > £ 0 ; 4> , j< ) 



-m CO — CO — OO — OO 



d^dKdrdr^ . 



(2.2-51) 



If we define the energy of the output signal 



E 



y 



as 



OO OO 




ao — OO 



| y ( t , r ) l^dtdr. 



(2.2-52) 



then the average output energy can be expressed as 



Ey & E<E y > = 



OO OO 

// 

— OO — OO 



E { |y(t,r) I }dtdr (2.2-53) 



or 



E = 

y 




R ( t , t , r , r ) dtdr , 



— OO OO 



(2.2-54) 



If Eq. (2.2-51) is substituted into Eq. (2.2-54), then 



OO OO OO 



E /E = 
Y x 




R„ ( x , r r i <) d^dKdxdr 

3) “O O 



— OO — OO — OO — . OO 



(2.2-55) 
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where 



E 

x 



CD OCX 




— OO — OO 



2 dtdr 



(2.2-56) 



is the energy of the input signal. Equation (2.2-55) indicates 
that the ratio of the output (received) average energy to the 
input (transmitted) energy for a WSSUS communication channel 
can be obtained by integrating the scattering function of the 
channel. Also note that the average output energy is not a 
function of the input signal's shape. 

Alternate expressions for the output autocorrela- 
tion function can be obtained from Eq. (2.1-25). If Eq. 

(2.1-25) is substituted into Eq. (2.2-47), then 



CD OO OO OO 



R y (t, t' ,r,r ' ) = 




X ( f , v ) X* (f ' ,v' ) 



— OO — OO — OO — OO 



Ry (f , f ' ,y_,_v ' ; t,t',r,r')* 

exp [ + j 2 tt (f t-f ' t ')] exp [- j2ir (yy r-y ' • r ' ) ] 

dfdvdf'dv'. (2.2-57) 



If the random filter is a WSSUS communication channel, then 
the autocorrelation function of the transfer function is given 
by Eq. (2.2-20). Substituting Eq. (2.2-20) into Eq. (2.2-57) 
yields 
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R < t, t' ,£,£') 



Co oo co ao 




X (f , v) X* tf ' /£' ) Rjj ( Af , A\^ ; At,Ar) 
exp[ + j2ir{ft - f ' t ' ) J exp [- j 2 tt ( \> *r 



-v* *r ') ] 



df dvdf ’ dv 1 . 



(2.2-58) 



A relationship between the input and output power spectral 
densities will be obtained next. 

Let us assume that the input signal x(t,r) is a 
zero-mean, wide-sense stationary (in time and space) , random 
process which is uncorrelated with the transfer function. 
Under these assumptions , the output autocorrelation function 
given by Eq. (2.2-57) becomes 



OO OO OO CO 



R (t,t‘ , r ,r ' ) 



- / / / / 



— OO -• OO « OO «QO 



E{X(f,v)X*(f ' ,v') }* 

Rjj (f / f ' / £/£' ; t , t ' /£/£') • 

exp [+ j 2n (f t-f 1 1 ')] exp [- j 2u ( £•£-£' •£' ) ] 

dfdvdf'dv' . (2.2-59) 



Since x(t,r) was assumed to be wide-sense stationary in both 
time and space, then it can be shown that 



E{X(f,v)X*(f', v') } = S (f , v) o (f-f 1 ) 6 (y-v' ) 

(2.2-60) 



where 



X (f , v) = F t F r (x(t,r) }, 
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S x (f,v) 4 F At F Ar {R x (At,Ar) } (2 . 2-61) 

is the power spectral density of the input, and 



R x (At,Ar) = E{x(t,£)x* (t' ,£’) } 



where At = t-t' and A£ = £-£' • If Eq. (2.2-60) is substituted 
into Eq. (2.2-59), then 



R y (t,t' ,£,£') 




^ x ( f , y.) R^j (f , f, v, v ; 



exp (-j 2ttv • Ar) dfdv , 



t, t',£, £')exp ( + j 2irf At) • 
(2.2-62) 



and if it is further assumed that the transfer function is wide- 
sense stationary in both time and space, then Eq. (2.2-62) re- 
duces to 



R y ( At, A£) 




S x (f,v) Rjj ( f , f , £, v_ 



At, A£)* 



exp ( + j 2-tff At) exp ( - j 2iry_* A£) df dy 



(2.2-63) 



which implies that the output y(t,£) is also wide-sense sta- 
tionary in both time and space. If we define the output power 
spectral density S y (n,B_) as 



S ( n , 8) 

y 






(2.2-64) 
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then substituting Eq. (.2.2-63) into the R.H.S. of Eq. (2.2-64) 
yields 



Sy (n # S.) 




( f , x j _) Rg (f,f,v,'J 



n-f, 8_-v_)dfdv_ 
(2.2-65) 



where, from Eq. (2.2-41), 



Rg(f,f,y_/v_ ; n-f , 3_-y_) 



oo oo 




— CO — oo 



R ( f , f , v , v ; A t, Ar) exp (-j 2:r (n-f ) At] * 

n — — — 

exp [ + j 2 7T ( S_-v_) • Ar ] dA tdAr . 



( 2 . 2 - 66 ) 



Note that Eq. (2.2-65) is in the form of a multidimensional 
convolution integral analogous to Eq. (2.1-46) for the deter- 
ministic case. The convolution process accounts for the fre- 
quency spreading (frequency in HZ. and spatial frequencies in 
cycles/m.) of the input power spectral density. 



Example 2.2-1 

Using Eqs. (2.2-65) and (2.2-66), let us compute the 
output power spectral density for the case when the transfer 
function is time-invariant, space- invariant, and deterministic . 
First note that 



R H (f,f,jj,v^ ; At,Ar) 



E { H { f , v ; t,r)H*(f,v ; t',r')}. 
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If H(jf,_v ; t,r) is both time- invariant and space- invariant , 
then [see Eq. (2.1-16)] 



H(f , v ; t,r) = H ( f , v) , 



and if the transfer function is also deterministic, then 



; At, Ar) = E{|H(f,v_)| } 



= | H ( f , v ) 



Therefore, Eq. (2.2-66) becomes 



Rg ( f , f , v_, ; n-f,S-jj) = | H ( f 



f-v.l 2 j 



exp [— j 2ir ( n-f ) At]dAt 



oo 

j exp [+ j 2 tt ( j3-v_) ♦ Ar ] 



dAr 



= | H ( f , v ) | 5 (n-f ) 6 (6-v), 



and substituting this expression into Eq. (2.2-65) yields the 
desired result 



S (n , 3) = | H ( n , 6 ) | S ( n , 3 ) . 

y — — a — 

By replacing n with f and 3. with v_, we finally obtain {compare 
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with Eq. (2.1-48)] 



S (f,v) = |H(f,v) | 2 S x (f,v) . 



k k k kk 
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III. COUPLING EQUATIONS 

We will now present those equations which couple the 
transmitted and received electrical signals to the transfer 
function of the ocean medium via the transmit and receive far- 
field directivity functions. Referring to Fig. 2, the fre- 
quency and angular spectra of the input acoustic signal to 
the medium and the output electrical signal from the receive 
aperture (array) are given, respectively, by the following 
equations [2] : 



where the frequency and angular spectrum of the output acoustic 
signal from the medium is given by [see Eq . (2.1-46)] 




(3-1) 



and 




oo oo 




x^(f , v_) B M (f ,y_ ; n-f , B_-y_) df dv 



(3-3) 
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where [see Eq. (2.1-32)] 



B M ^ f ' — ; 



GO 00 




; t,r) exp(-j2Tn>t) • 
exp ( + j 2ttk_* r) dtdr 



(3-4) 



is the bi-frequency function, and H M (f,y_ ; t,r) is the time- 
variant, space-variant, transfer function of the ocean medium. 
Equations (3-1) thru (3-4) are the basic coupling equations . 
However, w.r.t. Eq. (3-1), we have the following additional 
relationships : 



V f '^> = F t F r {x M (t '^ } 




x^(t,r) exp (+j 2 ttv_* r) dV« 
exp (- j 2uf t) dt (3-5) 



where [2] 






oo 

j X(f ,r) A T (f ,r) exp (+j2irft) df , 

— oo 



(3-6) 



X(f , a) = F t F r {x(t,r) } = 




x ( t, r) exp (+ j 2ira • r) dV- 



— 00 



exp ( - j 2irf t) dt 



(3-7) 
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is the frequency and angular spectrum of the transmitted elec- 
trical signal, and 



D T (f,a) = F r (A T (f,r)} = J A T (f ,r) exp ( + j 2ira *r) dV 

V 

(3-8) 

is the far-field directivity function (beam pattern) of the 
complex transmit aperture A T (f,r). Also, w.r.t. Eq. (3-6), 



oo 

X(f,r) = F t (x(t,r) } = j x ( t, r) exp (- j2irf t) dt 

— oo 

(3-9) 



and 



-1 



oo 

/ 



A T (f,r) = F ^ {D T (f ,a) } = j D T ( f , a ) exp (- j 2ira • r ) da . 

(3-10) 



Similarly, w.r.t. Eq. (3-2), we have the following 
additional relationships: 



Y(n,y) = 



F t F r {y } 



oo 

/ 



y ( t , r ) exp ( + j 2 tty/ r ) dV* 



exp (- j 2irnt) dt 



(3-11) 
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where [2] 



y (t,r) 



<30 

j Cn , r ) A R ( n , r) exp (+ j 2irnt) dn , 

— OO 



(3-12) 



the frequency and angular spectrum of the output acoustic sig 
nal from the medium is given by 



Y M (n,S) 



F t F r {y M (t '^ )} 




y (t,r) exp ( + j 2 tt S_* r) dV* 
exp (- j 2irnt) dt (3-13) 



where [see Eq. (2.1-25)] 



oo oo 



Y M (t ' 



-// 



— oo — oo 



(f ' — ) H m (f ' — ' t /£) exp ( + j 2nf t) 



exp (- j 2irv_- r) dfd\^ , 



(3-14) 



and 



D R (n,S_) = F r (A R (n,r)} = j A R ( n , r) exp ( + j 2nS_-r) dV 

V 

(3-15) 

is the far-field directivity function (beam pattern) of the 

complex receive aperture A (n,r). Finally, w.r.t. Eq. (3-12) 

R. — 



58 



co 

-r)} = J 



Y M ( n/£) = F t (y M (t,r)} = I y M Ct,r)exp(-j2un t)dt 



(3-16) 



and 



A R ( n ,r) = F" 1 {D R ( n ,g_)} = 



oo 

= / D 



R (n /B.) exp (- j 2 it§_* r) d§_. 

(3-17) 



Substituting Eq. (3-14) into Eq. (3-16) yields 



oo oo 



SI 

— oo —oo 



Y M (n,r) = I I X M (f ,v_)B M (f ,v_ ; n-f ,r) exp (- j2nv_- r) dfdv_ 



(3-18) 



where 



B M (f,v_ ; <p,r) = J H M (f,v_ ; t, r ) exp (- j 2 tt 4> t) dt 

— 00 

(3-19) 

which is not the same as the bi-frequency function [see Eq. 

(3-4) ] . 

The components of the vectors a, B_, and y_ are spatial 
frequencies with units of cycles/m. and r = (x,y,z) . 

A moments reflection leads one to the conclusion that 
Eq. (3-12) is well suited for space-time signal processing 
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applications. By inspecting Eqs, C3-18) , (.3-19) , and (3-1), 

it can be seen that the output signal given by Eq. (.3-12) can 
be expressed as a function of time and space in terms of the 
transmitted electrical signal, the transmit aperture, the 
transfer function of the ocean medium, and the receive aperture. 
Equation (3-2) , however, represents the theoretical output 
"spectrum", albeit dependent upon the same system functions. 

If the ocean medium is modelled as a random filter, 
then y(t,r) given by Eq. (3-12) is a random process with auto- 
correlation function 

Ry (t, t ' , r , r ' ) = E{y (t,r) y* (t' ,r ' ) } . (3-20) 

It is probably obvious by now that the weak link in the 
coupling equations is the ocean medium transfer function 
Hj^ ( f , y. ; t,r) and its corresponding autocorrelation function 
[see Eq. (2 . 2-2) ] 

Rjj (f,f',v,v' ; t,t',r,r’) = E{H M (f,v ; t,r)H*(f',v' ; t',r')}. 
M 

(3-21) 

If one is not able to specify a realistic functional form for 
H j/j(f ,— ' ^ ' t ^ ien t ^ ie equations presented in this section 

represent an exercise in the methods of linear system theory. 
This leads us to the main purpose of this paper which is to 
demonstrate the derivation of an ocean transfer function based 
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upon the W.K.B. approximation. The derivation is discussed in 
Section 4.1. 

Finally, note that the autocorrelation function given 
by Eq. (3-21) can be referred to as a generalized coherence 
function since it is a generalization of the two-frequency cor- 
relation function or two-frequency mutual coherence function 
based upon linear, time-variant filter theory as discussed by 
Ishimaru [51], for example. Equation (3-21) provides informa- 
tion concerning the amount of spreading in time-delay, space, 
frequency, and spatial frequencies that a transmitted signal 
will be expected to experience as it propagates in a random 
medium. The result of all this spreading is, of course, dis- 
tortion in pulse shape. 
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IV- ANALYSIS 

4 . 1 Transfer Function 

We will now proceed to derive a transfer function 
which models the bistatic communication channel geometry shown 
in Fig. 1. The communication channel is regarded to be the 
ocean volume between the apertures so that surface and bottom 
scattering effects are not included. Both apertures are sta- 
tionary (not in motion) , and it is assumed that no discrete 
point scatterers (such as bubbles, fish, etc.) are in the 
volume between the apertures. No motion implies that the re- 
sulting transfer function will be time- invariant. 

The propagation of small amplitude acoustic signals 
in the ocean from the transmit aperture to the receive aperture 
can be described by the following linear , inhomogeneous, scalar 
wave equation: 



V *(t,r) - 



c 2 (r) 



at 



2 *(t,r) = x M (t,r) 



(4.1-1) 



where is the velocity potential at time t and position 

r = (x,y,z), x..(t,r) is the source distribution [see Fig. 2 

— iU — 

and Eq. (3-6) ], and c (r) is the speed of sound in the ocean. 
Since the coupling equations discussed in Section III already 
allow for an arbitrary x M (t,r) with corresponding frequency 
and angular spectrum X M (f ,y ) , we need only find the solution 
to the following Helmholtz wave equation: 



62 



V 2 <?Cr) + k 2 n 2 (r)*>(.r) = 0 (4.1-2) 

where 

k = 2tt f / c = 2tt/a (4.1-3) 

is the constant, reference wavenumber, 

n(r) = c Q /c(r) (4.1-4) 

is the random index of refraction, 

c = c(r ) = fX (4.1-5) 

o—o o 

is the constant, reference speed of sound at the source posi- 
tion = (x o ,y o ,z o ) , and 

¥>(t,r) = <^(r) exp ( + j 2iTf t) (4.1-6) 

is the time-harmonic solution of Eq. (4.1-1) when x M (t,r) is 
set equal to zero, and where <P(r) is the solution of Eq. (4.1-2) 
Note that the wavenumber 

k (r) = 2irf/c(r) (4.1-7) 



can be expressed as 



k (r) = k Q n (r) . 



(4.1-3) 
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Therefore, k(r ) = k since n(r ) =1. The index of refraction 
— o o — o 

is commonly written as [40, 42-45J 



n(r) = n D (r) + n R (r ) (4.1-9) 

or 

n (r) = n D (r) + o(r)n NR (r) (4.1-10) 

where (r ) is the deterministic component and is usually 
close to unity in value, n R (r) is the random, zero mean compo- 
nent, air) is the standard deviation of n R (r) , and 

n^ R (r) = n R (r)/a(r) (4.1-11) 



is the normalized random component with zero mean and variance 
equal to unity. We shall work with Eq. (4.1-10) in this paper. 
Note that the average value of n (r ) is equal to n^ (r ) . 

Let us assume that the speed of sound is only a 
function of the depth y, i.e., c (r ) = c(y), so that Eq. (4.1-2) 
reduces to 



V <P(x,y,z) + 



k^ n' 2 (y) <?(x,y, z) = 0 (4.1-12) 



where, from Eq. (4.1-10), 



n(y) = n D (y) + a(y)n NR (y). 



(4.1-13) 
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By using the method of separation of variables and the W.K.B. 
approximation [46,47], an approximate solution of Eq. (4.1-12) 
is given by 



*>(x,y,z) s; exp[-jk x (x-x o ) ] a y (y ) exp [ + je y (v) ] exp[ - jk z ( z-z q ) ] 



where 



(4.1-14) 



a y (y) = Ek Y (y )]~ 1/2 



(4.1-15) 



e y (y) 



exp 



y 

“j j k y (?)dc 

y o 



(4.1-16) 



and 

where k^, k y (y), 
vector 



1/2 

k y (y) = t k o n2(y) " k x" k Z ] (4.1-17) 

and k z are the components of the propagation 



k(y) = k y x + k y (y)y + k z z. (4.1-18) 

Note that k y and k z are constants while k y (y) is a function of 
the depth y. In addition, Eqs . (4.1-14) thru (4.1-16) allow 

for a general source location with position vector = (x Q ,y o , z q ) 
The W.K.B. approximation given by Eqs. (4.1-15) and 
(4.1-16) is a valid solution for the depth dependence provided 



65 



that [46,47,50] 1) the transmitted frequencies f are high, 2) 
the sound speed profile c (y) is slowly changing, and 3) the 
depth interval from y Q to y does not include any turning 
points. A turning point exists at y = y T if k^(y T ) = 0. 

Since the reference propagation vector k can be 



expressed as 



then 



where 



JSq = k x x + k Y y + k z Z ' 



(4.1-19) 






+ k 



(4.1-20) 



k 



X 



k u 
o o 



(4.1-21) 



i 



and 



where 



k v = k v (4.1-22) 

Y O O 

k„ = k w (4.1-23) 

Z o o 

u = sine ccsij> (4.1-24) 

o o o 



and 



v = sine sinio = cosg (4.1-25) 

o o o o 

w = cose (4.1-26) 

o o 



are the direction cosines w.r.t. the positive X,Y, and Z axes, 
respectively, and (Qq/'I'q) are the vertical and azimuthal 
spherical angles measured w.r.t. the positive Z and X axes, 
respectively, representing the initial directions of wave pro- 
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pagation (see Fig. 7). From Eq. (4.1-20) we can write that 



, 2 . . 2 _ . 2 ,2 

k X + k Z k o ’ k Y 



(4.1-27) 



and substituting Eqs. (4.1-27) and (4.1-22) into Eq. (4.1-17) 
yields 



k Y (y> = ky 



1 + [n^(y)-11 



> 1/2 



(4.1-28) 



Note, that (y ) = k y since n(y Q ) = 1. Therefore, the square 
root expression in Eq. (4.1-28) is responsible for changing 
the initial direction of propagation k y . If 

| [n 2 (y ) - l]/v 2 | < 1 , (4.1-29) 

then the square root expression in Eq. (4.1-28) can be approx- 
imated by the first two terms in a binomial expansion yielding 

k Y (y) » k y + k 2 [n 2 (y) - l]/(2k y ) . (4.1-30) 

It will be shown later that the binomial expansion criterion 
given by Eq. (4.1-29) can be related to the critical angle of 
incidence, and hence, total reflection. The assumption con- 
cerning the absence of turning points will also be discussed 
later, but first let us return to the derivation of the trans- 



fer function. 
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Recall that we designated the vector v_ to represent 
the transmitted spatial frequencies as in X M (f,v_) (see Fig. 2) . 
If we let 

v = (f x ,f y ,f z ) , (4.1-31) 

then 

f v = k v /(2T7) = u/ X (4.1-32) 

A A O O 

f Y = k Y /(2TT) = v o /X o (4.1-33) 

and 

f z = k z / ( 2 tt ) = w o /X Q (4.1-34) 



are the spatial frequencies in the X,Y, and Z directions, re- 
spectively, where c Q = c(y Q ) = fX Q . Therefore, if Eqs. (4.1-30) 
thru (4.1-34) are substituted into Eqs. (4.1-14) thru (4.1-16), 
then Eq. (4.1-6) becomes 



exp 



-j [k Q / (47rf y ) ] 



/<■= 



(C)-I]d5> 



*(t,x,y, z) « 



j2 tt f Y + k 2 [n 2 (y) - 1] /(4irf Y )| 



1/2 



exp (+ j 2irf t) exp [ — j 2 tt (r-r^ ] . (4.1-35) 



Since the input to the communication channel is the time-har- 
monic plane wave 



exp ( + j 2irf t) exp [- j 2-n vy (r-r^ ] , 
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and since Eq. (4.1-35) represents the output from the channel 
at any time t and position r, then [see Eq. (2.1-21)] 



• t,r) = H M (f,v_ ; £> = H M ( f / f Y ; 

(4.1-36) 



where the random, time- invariant , space-variant transfer func- 
tion of the ocean medium is given by 



where 



H M (f ' f Y ; y) = ( f / f y ' y)exp[ + j0 M (f,f Y ; y) ] 

(4.1-37) 



V f ' f Y ; Y) = {2TTf y + k^[n 2 (y)-l]/(4Tf Y ) }" 1/2 

(4.1-38) 



^ 9 9 ^ 



Y 

-[k 2 /(4^f Y )] J [n 2 U)-l]d C 



(4.1-39) 



k = 2-irf/c 
o o 



(4.1-3) 



and from Eq. (4.1-13), 

n 2 (y) = n 2 (y) + 2n D (y ) a (y ) n NR (y ) + a 2 (y)n 2 R (y). 

(4.1-40) 



If the medium is "weakly irregular" or "weakly scattering". 
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i.e., if a (y) is very small compared to unity so that 



v ty) s Y << 1 



(4.1-41) 



where s^ is the scale size of the irregularities (i.e., the 

average distance over which the refractive index fluctuations 

2 

remain correlated [48]), then terms involving a (y) can be ne- 
glected [42] and Eq. (4.1-40) reduces to 



n 2 (y) a n 2 (y) + 2n D (y ) a (y) (y ) . (4.1-42) 



Therefore, by using Eq. (4.1-29) to simplify Eq. (4.1-38) and 
upon substituting Eq. (4.1-42) into Eq. (4.1-39), one obtains 
the following: 

Vf, f Y ; y) = ( 2 tt fy) -1/2 (4.1-43) 

and 

V f ' f Y ; y) = 0 MD (f ' f Y '* Y) + 0 MR (f ' f Y ; y) 

(4.1-44) 

where y 

0 MD (f ' f Y ; Y) = “[ k o /(4lTf Y^ J [ n D ( ^ )_ 1 J d? (4.1-45) 

y o 

is the deterministic or average component of the phase function. 



and 
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y 



0 MR Cf ' f Y ' = -[^/(27rf Y )] 



/ 

y o 



n D (c )a U ) n^U )dc 



(4.1-46) 



is the random component. 

The medium transfer function H given by Eq. (4.1-37) 
with amplitude and phase functions and 0 specified by Eq. 
(4.1-43) and Eqs. (4.1-44) thru (4.1-46), respectively, in- 
dicate that for a weakly irregular medium, the major effect of 
the medium is to angle modulate the transmitted field. Further- 
more, if n D (y) « 1, then 0 MD (f,fy ; y) « 0 [see Eq. (4.1-45)] 
and, as a result, the angle modulation is due strictly to the 
random fluctuations a(y)n^ R (y) of the index of refraction [see 
Eq. (4.1-46)]. Note that the angle modulation process is often 
referred to as "scattering" [48]. Also note that the transfer 
function derived in this section can be written as the product 
of two functions, one deterministic and the other random, i.e., 



which agrees with the assumed transfer function expressions of 
Laval and Labasque [36] and with the general practice of re- 
presenting a field in a random medium as the product of a de- 
terministic and a random function [40] . 



H. 



‘M 



[A M exp( + j0 JilD ) ] [exp( + j0 MR ) ] 



Let us now discuss the two important assumptions 



responsible for the derivation of the transfer function, namely, 
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the aforementioned binomial expansion criterion given by Eq. 
(4-. 1-29) and the assumption regarding the absence of turning 
points. 

Since 

n(y) = c Q /c(y), (4.1-47) 



substituting Eqs. (4.1-25) and (4.1-47) into Eq. (4.1-29) 
yields 



g < cos 
o 



-1 



(|c 2 (y) - c 2 |) 1/2 



c (y) 



(4.1-48) 



where g Q is the angle of incidence of the reference propagation 

vector (see Fig. 7) . Now recall that the critical angle of 

incidence 6 associated with a time-harmonic olane wave inci- 
c 

dent upon a plane boundary between two fluid media is given by 
[49] (see Fig. 8) 



sinB c = c i/ c 2 > c 2 > c i (4.1-49) 



or, equivalently. 



g = cos 



-1 



, 2 2 , 

(c 2 - c x ) 



1/2 



C 2 > C 1 



(4.1-50) 



where must be greater than c^ for B c to exist. When the 

angle of incidence g > g , there is total reflection and hence, 

c 
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Fig. 7. The reference propagation vector 

k and associated angles. 

— o 
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Fig. 8. Illustration of the angle of incidence 

8 and the critical angle of incidence B c » 
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no transmission of energy into medium II. The speed of sound 
c (y) generally increases with depth, for example, in the deep 
ocean. Also, c Cy 1 increases above and below the SOFAR channel 
axis. Therefore, when c(y) > c q , the absolute value sign in 
Eq. (4.1-48) can be removed, and by comparing Eq . (4.1-50) 

with the R.H.S. of Eq. (4.1-48), the binomial expansion crite- 
rion indicates that the initial angles of transmission 8 must 

o 

be less than the "critical angle" in order to avoid total re- 
flection, and thus, passing thru a turning point [47,50]. 

This is consistent with the fact that the W.K.B. approximation 
is invalid at a turning point [47,50]. 
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4 . 2 Output Electrical Signal 

Now that we have derived an ocean medium transfer 
function, let us demonstrate the use of the coupling equations 
presented in Section III by calculating the output electrical 
signal y(t,r) from the receive aperture (array) as given by 
Eq. (3-12) . As can be seen by inspecting the integrand of 
Eq. (3-12), we need expressions for both of the kernels Y,j(n,r) 
and A (ri/r). Let us compute Y (ri/r) first [see Eq. (3-18)]. 

Assume that the transmit aperture depicted in Fig.l 
is a planar array of MxN (odd) complex weighted point sources, 
centered at (x o ,y o ,z Q ) and parallel to the XY plane. In addi- 
tion, assume that the complex weights are separable. Since in 
most practical situations an identical input electrical signal 
is applied to all elements in the transmit array before the 
complex weights, i.e., since 

x ( t, r) = x ( t) , 

then Eq. (3-7) reduces to 

X(f ,o) = X(f) 5 (a) (4.2-1) 

since F^fl} = 6 (a_) . Substituting Eq. (4.2-1) into Eq. (3-1) 
yields 



X M (f,v) = X(f)D T (f,v) 



(4.2-2) 
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where 




exp (+j 2Trf y nd Y ) exp (+ j 2TTf x X Q ) - 
exp( + j2Trf Y y o )exp(+j2Trf z z o ) 



(4.2-3) 



is the far-field beam pattern of the transmit array, and d n 
are complex weights, d^ and d Y are the interelement spacings in 
the X and Y directions, respectively, and the last three expo- 
nentials are phase factors which account for the array being 
centered at (x o ,y Q > Z Q ) instead of at the origin (0,0,0). 

Since the transfer function derived in Section 4.1 
is time-invariant, Eq. (3-19) becomes 



B M (f,v_ ; <j> , r) = H M (f,v_ ; r)6(<f>) 



(4.2-4) 



and, as a result. 



B M (f,v_ ; l-f,r) = H M (f,v_ ; r)'S(n-f). 



(4.2-5) 



Therefore, substituting Eqs. (4.2-2) and (4.2-5) 



into Eq. (3-18) yields 
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/ 



Y M (n,r) = X(n) j D T (n ,v_) H M (n , v_ ; r) exp (-j2^v_* r) dv_. 



M 



(4.2-6) 



and upon substituting Eq. (4.2-6) into Eq. (3-12), we obtain 



Y (t,r) 



/ / 

— OO —00 



D T (f,v) H M (f,v 



r ) exp ( - j 2vv • r ) dv_* 



(f , r) exp ( + j 2 tt f t) df 



(4.2-7) 



where n was replaced by f since there is no frequency (HZ) 
spreading in this example due to the time- invariant property 
of the transfer function. We need to specify the complex 
aperture function A R (f,r) next. 

Assume that the receive aperture depicted in Fig. 1 
is a planar array of M'xN* (odd) complex weighted point 
sources, centered at (x R ,y R ,z R ) and parallel to the XY plane. 

In addition, assume that the complex weights are separable. 
Therefore, the receive aperture function is given by 



(M'-l)/2 (N 1 - 1 ) / 2 

V f '£) = S S c>d^o(x-[x R+ id']) • 

i=-(M’-l)/2 q=-(N’-l)/2 

My-[y R + qdy] ) 6 ( z-z R ) 



(4.2-8) 
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where c! and d' are complex weights and d' and d ' are the inter- 

1 q c X Y 

element spacings in the X and Y directions, respectively. 

Upon substituting Eqs . (4.1-31), (4.1-36), (4.2-3), 

and (4.2-8) into Eq, (4.2-7) and recalling that r = (x,y,z), 
one obtains 



(M'-l)/2 (N ' -1) / 2 

y ( t , x, y , z) = 2 c i d q 

i=-(M'-l)/2 q=- (N ' -1) / 2 

(M-D/2 

00 

2 c m f ex P'-J 2 ” f x iX im )df X- 

m=- (M-l)/2 -« 




(N-l)/2 

CO 

d n f V f ' f Y ; y R + q^)exp(-j27Tf Y AY qn )df Y 

n=- (N-l) / 2 
00 

J exp (- j 2trf Z AZ ) df ^exp (+ j 2nf t) df • 

— 00 

o(x-[x R + id^] ) 5 (y- [y R + qd^])6(z-z R ) (4.2-9) 

where 

AX im = (x R~ X o ) + (id X “ md X } (4.2-10) 

AY qn = (y R" y o ) + (qd Y ~ nc V (4.2-11) 

and 

AZ = z~ - z . (4.2-12) 

R o 
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If we now change variables from spatial frequencies to direc- 
tion cosines by substituting Eqs . (4.1-32) thru (4.1-34) into 

Eq. (4.2-9), and treat the frequency variable f as a constant 
with respect to the spatial frequency integrations, then Eq. 
(4.2-9) becomes 



(M'-l)/2 (N'-l)/2 



y (t,x,y, z) = (2/c^) 



i=-(M’-l)/2 q=-(N'-l)/2 






(M-l)/2 




c sine (2f AX. /c ) 
m lm o 



m=- (M-l)/2 



(N-l)/2 



+ 1 




n=- (N-l)/2 



a 



q 




(b exp(-jnb fAZ/c ) sine (b fAZ/c )}• 
q q O q O 



exp ( + j2^ft) df 6 (x- [x R + id£])S(y- [y R + qd^])5(z-z R ) 



(4.2-13) 



since in our example problem 
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-1 < u < +1 



o 



(4.2-14) 



a < v < +1 
q O 



(4.2-15) 



and 



0 < w < h 
o q 



(4.2-16) 



where 



a q = ( l n D (Y R + qd Y ] ’ 1 I )V2 



(4.2-17) 



b q = (1 - |n 2 (y R + qd^) - l|) 1/2 (4.2-18) 




(4.2-19) 



and n 2 ( • ) was replaced by n 2 ( • ) in Eqs. (4.2-17) and (4.2-18) 



The transfer function is given by Eqs. (4.1-37) and (4.1-43) 
thru (4.1-46) with the exception that the lower limit of in- 
tegration y q i n Eqs. (4.1-45) and (4.1-46) must be replaced by 
y Q + nd^ which is the Y coordinate of a point source in the 
transmit planar array. Note that if n^ ( • ) < 1, then 



Equation (4.2-13) is the desired result, i.e., it represents 
the random output electrical signal from each element in the 
receive array in terms of the transmitted electrical signal, 
the transmit and receive arrays, and the random transfer func- 
tion of the ocean medium. If y(t,x,v,z) given by Eq. (4.2-13) 



b q = n D (y R + ( J d Y ) * 



(4.2-20) 
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is complex, then simply take the real part to obtain the real 
output electrical signal. 
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4 . 3 Coherence Function 

Upon inspecting Eq. (4,2-13), it can be seen that 
the autocorrelation function of the output electrical signal 
will depend upon the autocorrelation function of the transfer 
function 



Rn (f,f',f Y ,f^ ; y,y ' ) = E{H M (f,f Y 
M 



y)H*(f',f£ ; y') } 
(4.3-1) 



which is also known as the coherence function. Substituting 
Eqs . (4.1-37), (4.1-39), and (4.1-43) into Eq. (4.3-1) yields 

Rjj (f,f',f Y ,f£ ; y,y') = ^2ir (f y f y ) 1//2 ] E | exp (+ j [KI (y ) +K ' I (y ' ) ] ) 

(4.3-2) 

where 

K = - k 2 /(4irf y ) (4.3-3) 

K' = +(k^) 2 / (4 uf y ) (4.3-4) 



y 

I(y) = J [n 2 U) - l]dc (4.3-5) 

*o 

y' 

[n 2 ( <;) - 1] d c 




Ky') 



(4.3-6) 
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k = 2 tt f /c 
o ' o 



(4.3-7) 



k' = 2irf'/c 
o o 



(4.3-8) 



2 

and n (•) is given by Eq. (4.1-42) where the expectation appear- 



ing in Eq. (4.3-2) is the characteristic function of the random 
quantity [KI (y) + K'I(y')]. If it is assumed that I (y) is a 
real Gaussian random process, which implies that the index of 
refraction is a real Gaussian random process, then Eq. (4.3-2) 
can be written as 





exp ( -E { 0^ R ( f , f y ; y) }/2) • 




exp ( -E { 0^ R ( f ' , f ; y')}/2) (4.3-9) 



where 



y Y 




n Q U ) n D ( c ' ) a ( e ) o ( ? ' ) . 



R U,£')d<;dc 
NR 



(4.3-10) 
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Y ' y' 



E{0 2 R (f\f^. ; y’)> = [Ck^) 2 /(2^)] 2 j j n D ( s) n D ( c 1 ) a U) 0 U ' ) • 



y y 

J o J o 



R (c,c')d 5 dc 
"NR 



(4.3-11) 



E { 0 MR (f , f Y ; y)9 MR (f , <4 ; y’)} = + [k 2 /(2Tif Y ) ] [(k^) 2 /(27rf^) ] 



y y ' 



J j n D tS)n I 



(5 , )o(c)o(c') 



y y 

■*0 -*o 



R (c/C')dcdc' 

Vr 



and 



(4.3-12) 



R n NR (y,Y ' ) = E{n NR (y)n NR (y ' } } ' (4.3-13) 



If it is further assumed that the deterministic component n D (y) 
of the index of refraction is equal to unity [40, 42-45], and 
that the random component n R (y) , and hence, the normalized 
random component n NR (y) is wide-sense stationary, i.e.. 



R (y,y‘) = R_ (Ay) (4.3-14) 

*NR n NR 



where Ay = y-y', then [see Eq. (4.1-45)] 



0 MD (f ' f Y '* Y) 



0 



(4.3-15) 
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and 



Q MD (f, ' f Y ? y ’ } = °' 



(4.3-16) 



and Eqs. (4.3-10) thru (4.3-12) become, respectively. 



2 -, 2 (y-y o ) 

:9 MR (f ' £ Y : y>! - <y-y 0 > [i^J / [ a - (7^r]\. s (5)as 



-(y-y 0 ) 



(4.3-17) 



2 (y’-y ) 



E * ®MR ( f ' 7 f y ; y,) ( = (y ' 



-v[Wl /"[*■ 



-(y'-y 0 ) 



and 



(y, - y o } 



R (5)dc 
NR 



(4.3-18) 



®MR ^ ' y ^ ®MR ^ '^y ' y ^ 



(k k* 0 ) 
o o 



(2irf y ) (2irf^) 



( y - y Q ) 



(y-y' ) 



/ 



R U)dc + (y • 

*NR 



-v / 



(y - y o ) 



(y-y' ) 



-(y'-y 0 ) 



R <5><lc + 
NR 



u 

/ 



CR (C)dC - 



n. 



NR 



-<y-y 0 ) 



(y- y o) 

(y-y') 



ZR (C)dC 

^R 



(4.3-19) 
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where a is the constant standard deviation of the Gaussian, 
zero mean, wide-sense stationary, random component n R (y) of 
the index of refraction. 
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V. SUMMARY and CONCLUSIONS 

A consistent notation, fundamental input-output relations 
and various time-space transformations for both deterministic 
and random linear, time-variant, space-variant filters have 
been established. The notation is consistent in the sense that 
all of the various input-output relations which are based upon 
the general theory will reduce to the classical relations of 
linear, time-invariant filter theory. These results should be 
of interest to persons involved in the general area of linear 
systems theory, and not only to those involved in underwater 
acoustics . 

With the use of the method of separation of variables 
and the W.K.B. approximation, a mathematical expression of a 
time-invariant, space-variant, random transfer function of the 
ocean medium was derived . The transfer function was time-in- 
variant instead of time-variant because motion was not consi- 
dered in the present derivation. The transfer function modelled 
the ocean volume between transmit and receive apertures (arrays) 
The ocean volume was characterized by a random index of refrac- 
tion (sound speed profile) which was a function of depth. The 
index of refraction was decomposed into deterministic and random 
components . 

In addition to the transfer function derivation, two 
example calculations were made. The first example demonstrated 
the use of the coupling equations and involved the derivation 
of a mathematical expression for the random output electrical 
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signal at each element in a receive planar array of complex 
weighted point sources. The output signals were expressed 
in terms of the frequency spectrum of the transmitted electrical 
signal, the transmit and receive arrays, and the previously 
derived transfer function of the ocean medium. The first exam- 
ple demonstrated that an output electrical signal could be 
derived in a logical and straightforward fashion. The second 
example involved the derivation of the coherence function, i.e., 
the autocorrelation function of the transfer function. In 
order to obtain somewhat simplified results, it was necessary 
to assume that the random component of the index of refraction 
was Gaussian and wide-sense stationary. 
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